Phase space gaps and ergodicity breaking in systems with long range interactions 
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We study a generalized isotropic XY-model which includes both two-spin and four-spin mean- 
field interactions. This model can be solved in the microcanonical ensemble. It is shown that in 
certain parameter regions the model exhibits gaps in the magnetization at fixed energy, resulting in 
ergodicity breaking. This phenomenon has previously been reported in anisotropic and discrete spin 
models. The entropy of the model is calculated and the microcanonical phase diagram is derived, 
showing the existence of first order phase transitions from the ferromagnetic to a paramagnetic 
disordered phase. It is found that ergodicity breaking takes place both in the ferromagnetic and 
the paramagnetic phases. As a consequence, the system can exhibit a stable ferromagnetic phase 
within the paramagnetic region, and conversely a disordered phase within the magnetically ordered 
region. 

PACS numbers: 05.70.Fh Phase transitions: general studies 



I. INTRODUCTION 



The aim of this paper is to shed light on a dynamical 
feature which distinguishes systems with long range in- 
teractions from those with short range ones. It is well 
known that the attainable region in the space of exten- 
sive thermodynamic parameters is always convex when 
only short range interactions are present. Consider, for 
example, two magnetic subsystems with the same poten- 
tial energy V but with two different magnetization values 
mi and mi (notice that, instead of the magnetization, 
we might have chosen any other order parameter). In- 
troducing a parameter A which depends on the relative 
size of the subsystems, it is possible to get any interme- 
diate value of the magnetization m by just combining 
the two subsystems with appropriate weights, such that 
in = \mi + (1 — A)m2, while the potential energy is kept 
constant to V. It is important to stress that this con- 
vexity property is satisfied only if the system is additive, 
since the interaction energy between the two subsystems 
has been neglected. This can safely be done only for 
short range interacting systems in the large volume limit. 
Moreover, convexity implies that the accessible region in 
the space of thermodynamic parameters is connected. 

In contrast, systems with long range interactions [l[ are 
not additive. As a consequence, the convexity property 
may be violated, and then also connectivity. This fea- 
ture has profound consequences on the dynamics. Gaps 
may open up in the order parameters space. Such gaps 
have been recently reported in a class of anisotropic XY 
models [U, H| and for a discrete spin system Q ■ Since the 
accessible region in the order parameters space is no more 
connected, ergodicity breaking naturally appears when a 
continuous microcanonical dynamics is considered. This 
is the property we discuss in this paper. More specif- 



ically, the question we would like to address is that of 
locating the region in the phase diagram of the model 
where we expect ergodicity breaking to occur. We will 
argue that ergodicity breaking occurs in regions of pa- 
rameters where metastable states exist. In such regions, 
one may also generically expect first order phase transi- 
tions. Thus, we conclude that the breakdown of order 
parameters space connectivity is generically associated 
with first order phase transitions. Moreover, we con- 
jecture that the intersection between first order phase 
transition lines and the boundary of the parameter re- 
gion where the order parameters space is disconnected, 
occurs at non smooth (singular) points of the boundary 
itself. We also conjecture that the same is true for the 
mctastability line. 

In recent years, models describing classical rotators 
with all-to-all mean-field interaction have been inten- 
sively studied in order to explain general dynamical fea- 
tures of long-range interacting systems The most 
celebrated of such models is the so-called Hamiltonian 
Mean Field (HMF) model, where rotators denoted by the 
angle 6^ , i = 1 , . . . , N arc coupled through a long-range 
mean-field potential 
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The interaction is ferromagnetic for positive J-values, 
and antiferromagnetic otherwise. Defining the complex 
order parameter (with i the imaginary unit) 



m = (m, + im y ) 



TV ^ 



exp(i6>„), 



(2) 



the potential per particle V = W/N can be rewritten as 



-Jm 2 /2 where 



is the modulus of the magneti- 



zation. The total energy is finally obtained by adding a 
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kinetic energy term to the above potential so that we get 
the Mean-Field Hamiltonian 



N 2 

=1 



H = > f£ + JV V. 



(3) 



In this formulation, the model can also be thought of as 
describing a system of particles with unit mass, interact- 
ing through the mean-field coupling (the names rotator 
and particle might cquivalently be used). Note that a 
system with non-unit mass may be reduced to the above 
system by rescaling time. As the parameter J might be 
removed by a suitable energy rescaling, we will omit it 
in the remainder of the paper. Within the Hamiltonian 
framework, both the energy density E = H/N and the 
total momentum P = ^2 n p n are conserved quantities 
and they are fixed by the initial conditions. 

In this paper, we study a generalized mean-field Hamil- 
tonian Q with potential 
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which reduces to the antiferromagnetic HMF model (J = 
— 1) in the particular case A = 0. In addition to the 
first antiferromagnetic term, we consider a ferromagnetic 
fourth order term in the magnetization, whose intensity 
is controlled by a positive constant A. In the regime 
of large values of the parameter A, a magnetized state 
would be favorable because of the ferromagnetic coupling 
while, for small values of A, the antiferromagnetic cou- 
pling dominates and leads to a non-magnetic state. As 
we shall demonstrate below, there exists a range of the 
parameter A for which the model exhibits a first order 
phase transition between a paramagnetic phase at high 
energies and a ferromagnetic phase at low energies. In 
both phases there are regions in the (A, A) plane in which 
the accessible magnetization interval exhibits a gap, re- 
sulting in breaking of crgodicity. 

The organization of the paper is as follows. In Sec- 
tion HH we first address the issue of the accessible mag- 
netization interval in this model. In the following Sec- 
tion IIIIl we present the statistical mechanics of the 
model. In particular, we study the phase diagram and 
the conditions for phase transitions to occur. In Sec- 
tion IIVI we study the dynamical consequences of this 
crgodicity breaking. Finally, Section [V] will be devoted 
to conclusions and perspectives. 



II. PHASE SPACE CONNECTIVITY 

Let us study the phase space structure and the phase 
diagram of this model. The specific kinetic energy Ek = 
E—V (to) is by definition a non negative quantity, which 
implies that 



E > V{m) = m 2 /2 - Ato 4 /4. 



(5) 



We will show that as a result of this condition not all the 
values of the magnetization m are attainable in a certain 
region in the (A, A) plane; a disconnected magnetization 
domain might be indeed a typical case. As explained be- 
low, this situation is the one of interest. Let us charac- 
terize the accessible domains in the (A, A) plane more 
precisely by analyzing the different regions. 

For A < 1, the local maximum of the potential en- 
ergy V is not located inside the magnetization interval 
[0, 1] (see Fig. [IJi). The potential being a strictly in- 
creasing function of the magnetization, the maximum is 
reached at the extremum m = 1. The complete interval 
[0, 1] is thus accessible, provided the energy A is larger 
than V(l): the corresponding domain is in Rl defined 
in Fig. [5J The horizontally shaded region is forbidden, 
since the energy is lower than the minimum of the po- 
tential energy V(0) = 0. Finally, the intermediate region 
0<A<y(l)in included in the region R2: it is impor- 
tant to emphasize that only the interval [0, to_ (A, A)], 
where m± (A, A) = [(l± x /l-4AA)/A] 1 /2 ) 

is accessi- 
ble. Larger magnetization values correspond to a poten- 
tial energy V(m) larger than the total energy A, which 
is impossible. Figure [TJi also displays the value to_ cor- 
responding to an energy in the intermediate region R2. 

For A > 1, the specific potential energy V has a max- 
imum Vmax = 1/AK which is reached at < TO max = 
1/VK < 1. Figures [TJd and [TJ:, where the potential en- 
ergy ((4]) is plotted vs the magnetization m, present such 
cases. For an energy A larger than the critical value V mElfX , 
condition ([5]) is satisfied for any value of the magneti- 
zation to. The complete interval [0, 1] is thus accessible 
for the magnetization to. This region is Rl represented 
in Fig. H 

Let us now consider the cases for which A < V max . 
As discussed above, the minimum V mm of the potential 
energy is also important to distinguish between the dif- 
ferent regions. For 1 < K < 2 (sec Fig. [I}}), the minimum 
of V(m) corresponds to the non-magnetic phase m = 
where 1^(0) = 0. The quadrillcd region shown in Fig. [2j 
which corresponds to negative energy values, is thus non- 
accessible. On the contrary, positive energy values are 
possible and correspond to very interesting cases, since 
only sub-intervals of the complete magnetization inter- 
val [0,1] are accessible. There are however two different 
cases: 

• for < V(l) < A < Vmax, the domain of possible 
magnetization is [0, to_ (A, A)] U [m+ (A, A) , 1]. 
The above conditions are satisfied in R3 of Fig. [5] 

• for < A < V(l), only the interval [0,m_ (A, A)] 
satisfies condition ((5|). This takes place within R2 

of Fig. m 

For the domain 2 < A, the minimum of the potential 
energy is attained at the extremum, to = 1, implying 
A > V(l) = 1/2- A/4. The vertically shaded region is 
thus forbidden. In the accessible region two cases can be 
identified: 
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FIG. 2: The (E, K) plane is divided in several regions. The 
solid curve corresponds to K — 1/(4E), the oblique dashed 
line to K = 2 - 4E, while the dotted one to K = 1. The 
vertically shaded, quadrilled, and horizontally shaded regions 
are forbidden. The accessible magnetization interval in each 
of the four regions is indicated (see text for details) . 
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FIG. 1: Specific potential energy V vs. magnetization m for 
three different cases: K < 1 (panel a), 1 < K < 2 (panel b) 
and 2 < K (panel c). The location of the maximal magneti- 
zation mmax and the corresponding potential energy Vmax are 
shown (see text). In panels (a) and (c), two examples of the 
location of the critical magnetization m±(E,K) is indicated 
for energy values E in the intermediate regions. 



• for V(l) < E < 0, only the interval [m+ (E, K) , 1] 
satisfies condition ([5]). It is important to note that 
m+ (E, K) < 1 provided E > 1/2 - K/4. These 
cases correspond to R4- 



• for < E < Vinax, the two intervals [0, m- (E, K)] 
and [m+ (E, K) , 1] satisfy condition ([5]), corre- 
sponding to R3 of Fig. 



In summary, the complete magnetization interval [0, 1] 
is accessible only in the regioni?7. In R2, only [0, m_] 
is accessible, while only [m+, 1] is accessible in R4- Fi- 
nally, we note that the phase space of the system is not 
connected in the region R3. Indeed, the magnetization 
cannot vary continuously from the first interval [0,m_] 
to the second one [m+, 1], although both are accessible. 
These restrictions yield the accessible magnetization do- 
main shown in Fig. [3l The fact that for a given energy 
the phase space is disconnected implies ergodicity break- 
ing for the Hamiltonian dynamics. It is important to 
emphasize that the discussion above is independent of 
the number of particles and ergodicity is expected to be 
broken even for a finite (but sufficiently large) N. 
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FIG. 3: Accessible region in the (m, E) plane for K = 3. For 
energies in a certain range, a gap in the accessible magneti- 
zation values is present and denned by the two boundaries 
m±{E,K). 



III. STATISTICAL MECHANICS 

We have thus found that in certain regions in the 
(E, K) plane, the magnetization cannot assume any value 
in the interval [—1, 1]. For a given energy there exists a 
gap in this interval to which no microscopic configuration 
can be associated. In this Section, we study the statis- 
tical mechanics of model by considering the micro- 
canonical measure 



JY 



H{E, N) = \{ A6 k e\ Pk 5 (H - NE) 



(G) 



fe=i 



The probability P(m) to observe the magnetization m 
can be obtained for large values of the number N of par- 
ticles using large deviation techniques, in a similar way 
as it was derived for the HMF model [||. Here again, 
the appropriate global variables are the magnetization m 
and the mean kinetic energy Ek = J2nPn/N- Cramer's 
theorem allows one to derive the probability distribution 
P(m) in the large iV-limit. The entropy per particle, for 
given energy and magnetization, is obtained as 



s(E, to) 



lim — In P(m), 



(I) 



Following the steps described in Ref. [||, one finds that 
the specific entropy s can be expressed as 



s (E, to) = s K (E, to) + s con f (to) , 



(8) 



i.e. the sum of the momentum entropy, which is related 
to the kinetic energy 

s K (E, to) = In (E c ) 12 = \ In (e - ^ + , (9) 



and of the configuration entropy which, as for the HMF 
model, is 

Sconf (m) = -m<f>(m) + In (I (0(to))) . (10) 

Here <p is the inverse of the function I\ / Iq where I n is the 
modified Bessel function of order n. The microcanoni- 
cal thermodynamics is finally recovered solving the vari- 
ational problem 



S(E) = sup [s(E,m) 



(11) 



First, by comparing the low and high energy regimes, 
it is possible to show that a phase transition is present 
between the two regimes. In the domain R4 of Fig. [21 for 
very low energy E (close to the limiting value 1/2 — A"/4), 
the accessible range for to is a small interval located close 
to to = 1 (see Fig. [TJ;). The maximum of the entropy s 
corresponds therefore to a magnetized state, m eq , located 
very close to to = 1 (see the top- left inset in Fig. HJ. On 
the contrary, in the very large energy domain, the vari- 
ations of entropy with respect to to are dominated by 
the variations of the configurational entropy since the 
kinetic entropy (see Eq. © with to of order one and 
a very large energy E) is roughly a constant around 
s(E, 0) = {hiE)/2. As expected, the configurational en- 
tropy is a decreasing function of the magnetization: the 
number of microstates corresponding to a non-magnetic 
macrostate being much larger that the same number for a 
magnetized state. The configuration entropy has there- 
fore a single maximum located at m eq = 0. A phase 
transition takes place between the non-magnetic state 
at large energy, and a magnetic state at small energy. 
Moreover, as the non-magnetic state is possible only for 
positive energies E (see Figs.QJ, the transition line is lo- 
cated in the domain E > 0. In this region, the quantity 
d^ n s(E, 0) = —(1 + l/(2E)) is negative, which ensures 
that, for any value of E and K, the non-magnetic state 
to = is a local entropy maximum. The latter argu- 
ment allows to exclude a second order phase transition 
at a positive critical energy, since the second derivative 
d^ n s{E, 0) would have to vanish, which is impossible. The 
above argument leads to the conclusion that the phase 
transition must be first order. 

Let us now focus on the behavior of the entropy in the 
region R3, where the accessible range for to is the union of 
two disconnected intervals [0, to_] U [to + , 1]. As discussed 
above, the total entropy s(E, to) has a local maximum in 
the first interval [0, to_] located at to = and associated 
with the entropy = s(E : 0) = log(i?)/2. In the 

second interval [to + , 1] , a maximum is also present with 
s max = s(E,m eq ) where m eq > m + > 0. As s 1 l nax (i?) 
diverges to — oo when E tends to 0, a magnetized state 
is expected on the line E = 0, as long as s^ ax remains 
finite. Since K = 2 is the only value for which s^ ax (0, K) 
diverges, the first order transition line originates at the 
point B(0, 2) in Fig. 2J Although it is possible to study 
analytically the asymptotic behavior of the transition line 
near this point, we can rather easily compute numerically 
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the location of the first order transition line, represented 
by the dash-dotted line in Fig. 








-0.25 0.00 E 0.25 0.50 

FIG. 4: Phase diagram of the mean field model Q. The 
dash-dotted curve corresponds to the first order phase tran- 
sition line, issued from the point B(0,2). As in Fig. [2] the 
solid curve indicates the right border of the region R3, where 
the phase space is disconnected. The dashed line corresponds 
to the K = 2 — AE. The dotted line issued from the point 
.A(l/4, 1) represents the metastability line for the magnetized 
state, while the E = vertical dotted line is also the metasta- 
bility line for the non-magnetic state. The four insets repre- 
sent the entropy s versus the magnetization m for the four 
energies: E =-0.05, 0.005, 0.1 and 0.35, when K = 3. 

Figure [J] also shows the metastability line (the dot- 
ted line starting at point .4(1/4,1)), for the magnetized 
state m eq ^ 0. To the right of this metastability line, 
there is no metastable state (local entropy maximum for 
any m > 0, see bottom-right inset in Fig. [J} while a 
metastable state (local maximum) exists at some non 
vanishing magnetization on the other side (see top-right 
inset in Fig. HJ). Finally, the vertical dotted line E = 
corresponds to the metastability line of the non-magnetic 
state m = 0. 

One of the key issues we would like to address is the 
possible links between the breakdown of phase space con- 
nectivity and thus ergodicity breaking, on the one hand, 
and the phase transition, on the other hand. Obvious 
general properties do exist: a region of parameters where 
the phase space is disconnected corresponds to a region 
where metastable states do exist. Let us justify this state- 
ment. At the boundary of any connected domain, when 
the order parameter is close to its boundary value nib, 
there is a single accessible state. In a model with con- 
tinuous variables, like the one we study in this paper, 
this leads to a divergence of the entropy. In this case the 
singularity of the entropy is proportional to ln(m — m&) 
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FIG. 5: Entropy s versus magnetization m for K — 3 and 
several energy values. The different curves correspond, from 
top to bottom, to E = 0.35 (dotted), 0.155 (dash-triple 
dotted, metastability limit for the magnetized state), 0.1 
(dashed), E = 1/(4A") = 1/12 (solid, appearence of the gap), 
0.0089 (dash-dotted, first order phase transition), -0.05 (long 
dashed). This picture demonstrates that gaps in the accessi- 
ble states develop as the energy is lowered. 



(see for instance equation ©). For a model with discrete 
variables, like an Ising model, the entropy would no more 
reach — oo as m tends to mj, but would rather take a fi- 
nite value. However, the singularity would still exist and 
would then be proportional to (m — nib) hi(m — nib), hi 
both cases, of discrete and continuous variables, at the 
boundary of any connected domain, the derivative of the 
entropy as a function of the order parameter tends gener- 
ally to ±oo. As a consequence, entropy extrema cannot 
be located at the boundary. Thus, a local entropy max- 
imum (metastable or stable) does exist in a region of 
parameters where the phase space is disconnected. 

Hence, there is an entropy maximum (either local or 
global) in any connected domain of the phase space. For 
instance, considering the present model, in Fig. 31 the 
area R3 is included in the area where metastable states 
exist (bounded by the two dotted lines and the dashed 
line). In such areas where metastable states exist, one 
generically expects first order phase transitions. Thus the 
breakdown of phase space connectivity will be generically 
associated to first order phase transition, as exemplified 
by the present study. However, this is not necessary, 
one may observe metastable states without first order 
phase transitions, or first order phase transitions without 
connectivity breaking. 

A very interesting question is related to the critical 
points A and B shown in Fig. 0J As observed in the 
phase diagram, the end point for the line of first order 
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phase transition (point B) corresponds also to a point 
where the boundary of the region where the phase space 
is disconnected is not smooth. Similarly, the end point for 
the line of appearance of metastable states (point A) is 
also a singular point for the boundary of the area where 
the phase space is disconnected. It is thus possible to 
propose the conjecture that such a relation is generic, and 
that it should be observed in other systems where both 
first order phase transitions and phase space ergodicity 
breaking do occur. 



IV. DYNAMICS 

The feature of disconnected accessible magnetization 
intervals, which is typical of systems with long-range in- 
teractions, has profound implications on the dynamics. 
In particular, starting from an initial condition which 
lies within one of these intervals, local dynamics is un- 
able to move the system to a different accessible interval. 
A macrostate change would be necessary to carry the sys- 
tem from one interval to the other. Thus the ergodicity 
is broken in the microcanonical dynamics even at finite 
JV. 

In Ref. [i[ this point has been demonstrated using 
the microcanonical Monte-Carlo dynamics suggested by 
Creutz Q . Here, we use the Hamiltonian dynamics given 
by the equation of motions 

dH „ T / dm odm\ , . 

= TV (l — Km 2 ) (sin 6 n m x — cos 6 n m y ) . (13) 

We display in Fig. [6] the evolution of the magnetization 
for two cases, since we have shown above that the gap 
opens up when E decreases. The first case corresponds 
to the domain Rl, in which the accessible magnetization 
domain is the full interval [0, 1]. Fig.UJi presents the time 
evolution of the modulus m of the order parameter (J5|) 
versus time. The magnetization switches between the 
metastable state m = and the stable one m eq > 0. 
This is possible because the number of particles is small 
(JV = 20) and, as a consequence, the entropy barrier (see 
the inset) can be overcome. Considering a system with a 
small number of particles allows to observe flips between 
local maxima, while such flips would be less frequent for 
larger JV values. 

In the other case, we consider a stable m = state 
which is disconnected from the metastable one. It makes 
the system unable to switch from one state to the other. 
Note that this feature is characteristic of the microcanon- 
ical dynamics, since an algorithm reproducing the canon- 
ical dynamics would allow the crossing of the forbidden 
region (by moving to higher energy states, which is im- 
possible in the microcanonical ensemble). The result of 
two different numerical simulations is reported in Fig.j6j3. 
One is initialized with a magnetization within [0, m_], 
while the other corresponds to an initial magnetization 



a) 

1.0 F 




5000 
t 

FIG. 6: Time evolution of the magnetization m (the entropy 
of the corresponding cases is plotted as an inset). Panel a) 
corresponds to the case E — 0.1 and K = 8, while panel b) to 
E = 0.0177 and K = 3. In panel b), two different initial con- 
ditions are plotted simultaneously: the solid line corresponds 
to m(t = 0) = 0.1 while the dashed line to m(t = 0) = 0.98. 
The dashed (resp. dotted) line in panel b) corresponds to the 
line m = m+ ~ 0.794 (resp. m = m_ ~ 0.192). 



close to to = 1 (i.e. within [m+, 1]). One clearly sees 
that the dynamics is blocked in one of the two possible 
regions, and not a single jump is visible over a long time 
span. This is a clear evidence of ergodicity breaking. 



V. CONCLUSIONS 

We have found that for sufficiently low energy, gaps 
open up in the magnetization interval, to which no mi- 
croscopic configurations corresponds. Thus the phase 
space breaks into disconnected regions. Within the mi- 
crocanonical dynamics the system is trapped in one of 
these regions, leading to a breakdown of ergodicity even 
in finite systems. 

Ergodicity breaking implies that an ensemble average 
will not be well reproduced by a time average, usually 
simpler to get from the experimental or numerical point 
of view. Such a discrepancy, forbidden a priori for short- 
range interacting system, might be a typical case for sys- 
tems with long-range interactions. 

Of high interest would be of course to study systems 
with both short and long-range interactions since they 
usually compete in standard condensed matter systems 
@. Such a study could suggest an experimental system 
where ergodicity breaking might be observed. 
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